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Abstract 

In this paper is presented the infinitesimally locally Lorentzian theory of gravitation (ILLTG) , 
which predicts the same experimental results which confirm the GR. Motion in a neighborhood of 
the center of gravitation is considered. The geodesies between the ILLTG and GR are compared. 
ILLTG differs from GR mainly in the equations of motion and for v « c the trajectories are 
almost the same in both theories. Two tests in the solar system are proposed, such that ILLTG 
and GR predict different results. The mass and energy with respect to ILLTG is considered. The 
periastron shift of the binary stars is calculated and a cosmological problem is considered. 



1 Introduction 

In the paper [8] was presented a new model of gravitational theory, which yields to the same ex- 
perimental results for motion of a particle in a gravitational field predicted by the general relativity 
(GR) and which are experimentally verified. In this paper that theory is studied now in many details 
and a lot new results are presented. In this section will be repeated the basic concepts of the theory 
given in [8] with some improvements. We convenient to use that x'^ — ict for the time coordinate and 
orthogonal matrix A with complex elements will mean that AA^ = I. 

First we underline the main difference between the GR and the proposed theory. Let us start from 
the metric in GR 

ds^ = giidX"^ + g22dY^ + g^adZ^ + giid{icTY, 
where gij can approximately be written in the form 

V V 

gil = ff22 = .g33 = 1 + 2^, 544 = 1-2-2, gij = [i^i). 

Here V is the gravitational potential, for example V = in case of gravitation of a spherical body. 
Note that the coefficients g^ may not be determined very precisely like in the GR, because we shall 
use them for metric proportions, but not for equations of motion. Thus, the expressions of type c""* 
and less are neglected. Hence, the experiments for the red shift of the spectral lines and the Shapiro 
time delay, where the Christoffel symbols are not used, are explained in the same way like in the GR. 
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The above coordinates X, Y, Z, T will be called ordinary coordinates. In order to consider the 
equations of motion, we introduce the so called orthonormal coordinates x, y, z, t such that 



dx 



J—dX^ll-X^^dX, dy = x —dY^(l-X^)dY, rf^ = « fl - ^"jdZ, 

V 511 ^ V 522 V c^^ V 533 ^ 



Further we shall deal mainly with orthonormal coordinates, otherwise it will be emphasized that 
the coordinates are ordinary. According to the ordinary coordinates, the velocity of the light is 
equal to c^l — 2^^, like in GR, while according to the orthonormal coordinates the velocity of the 

light is universal constant. Indeed, according to the orthonormal coordinates the space-time in an 
infinitesimally small neighborhood is flat like in the special theory of relativity (STR), and hence the 
proposed theory will be called infinitesimally locally Lorentzian theory of gravitation (ILLTG). Note 
that ILLTG, using orthonormal coordinates, does not describes only gravitation, but also inertial 
forces (see Ref. [9]). The ordinary coordinates X,Y,Z,T from ILLTG are analogous to the usual 
coordinates x, y, z, t from the GR, while the orthonormal coordinates from ILLTG mean that the 
observer is far from gravitation, where the metric is like in the STR. 

Instead of the metric gij we can use only a potential ji. For example, this potential caused by a 
spherical body with mass M, whose center is at the coordinate origin and which rests with respect to 
the chosen coordinate system is /it = 1 + Thus, the equations (1.1) will be replaced by 

dx = iJ.~^dX, dy = fi'^dY, dz = fi'^dZ, dt = fidT. (1.1') 

Further we shall present the basic results about the equations of motion in orthonormal coordinates, 

and naturally then we should return back in ordinary coordinates as natural. This transformation in an 
infinitesimally small neighborhood is always possible, but in large we can deal with the nonholonomic 
coordinates. 

If two gravitational fields are given by the gravitational potentials /ii and /i2, then the superposition 
of those two fields is a gravitational field with potential /Xi/i2- In special case if IJ.2 = C is a constant 
in the space-time, then the motions of particles should be the same like for the field given by fii. 

ILLTG is built using a 4- vector potential Ui. If the magnitude of this vector changes but its 
direction is a constant, then there exists a gravitational field and the source of gravitational field is 
moving with 4-vector of velocity ?7i/||?7||; if the direction of this vector changes but its magnitude is 
a constant, then there exist some inertial forces, for example centrifugal. This 4-vector field can be 
written in the form 

(m,U,,Us,U,) = ^=^(^,^,^,l), (1.2) 

y^l - U^/C^ V «C IC IC I 

where u = (wi, U2, M3) is a 3- vector of velocity and \x is the gravitational potential. 

Example 1. The STR is characterized by /x = const, and u is a constant vector field. 

Exam,ple 2. The gravitational field of a spherical body with mass M which is at the origin of the 
coordinate system is given by the vector field (0, 0, 0, /x), where /U = 1 + 

Example 3. If we measure with respect to a noninertial system which rotates around the x^ axis 
with angular velocity w, then the 4-vector potential obtains by putting = 1 and u = {x^iv, —x^ui, 0) 
[9]. 

Further we introduce a non-linear connection induced by the 4-vector potential Ui. Let us denote 

by Sjds the form of connection in the direction of the unit vector V = dx'^/ds. The necessary and 
sufficient condition for the metric tensor gij = 6ij to be parallel, is 

Sj = -Si. (1.3) 
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Now let us introduce the following tensor 

1 



^{dUi/dx^ -dUj/dx'). 
If we denote by Vj the 4-vector of velocity of the moving particle, i.e. 

V2 V3 



^IC\/1 — v-^ 



one can easily verify that the following simple equations 

dV^jds = F.jVj, 



(1.4) 



(1.5) 



(1.6) 



which are analogous to the Lorentz force in the theory of electromagnetism, contain inertial accelera- 
tions (Coriolis acceleration, centrifugal and transversely acceleration) and gravitation acceleration for 
Newton approach. In order to obtain a precise theory, we consider the following two tensors 



1 



P-=S- 



fj, -]- U s y s h*' h*' 



(1.7) 
(1.8) 



These two tensors together with Fij are invariant under the gauge transformation ji C/i, where C 
is a constant. The tensors Fij and (pij are antisymmetric and Pij is an orthogonal matrix. It can be 
verified by using the identities C/ji7, = jj^ and ViVi = 1. Now the required connection Sij is introduced 

by 

Sij = Pri^rk-^kj ^ (-^''^) 

or in matrix form 5* = —P*(j)P. Although this connection seems to be very complicated, now we will 
show the opposite. Indeed, assume that (J7,) = (0,0,0, 1) at the considered point. Then 



F = 



where 



—liuj-^jc 2iijj2/c — ai/c^"' 

Hijj^jc — 2ia;i/c —a2/(? 

—2iw2/c 2iuji/c —aj,/c^ 

ai/(? a2/(? as/(? 

—iLOzjc iijJ2/c —ai/c?' 

iws/c —iuJi/c ~a2l<? 

—iui2/c iuii/c —azjc? 

a\l(? a2/(? azl<? 



a = (ai, a2, as) = c^grad + d\i/dt and w = —-^^oi u 



(1.10) 



(1.11) 



(1.12) 



represent the 3- vector of acceleration and the 3- vector of angular velocity. In that special case {JJi) = 
(0, 0, 0, 1), the tensor P is given by 



(1.13) 





-lVrV2 






-IV2V, 


1 - \yi 


-\V2Vi 


V2 




-hV3V2 


1 - 




. -Vi 


-V2 




V4. 
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where Vi, V2, V3, 14 are given by (1.5) and v = I + V4. Note that the matrix P represents just a Lorentz 
transformation for the 3-vector — v = {—Vi,—V2,—V3). Now the connection Sij becomes much more 
clear and acceptable. 

Remark. Although is a tensor, its components have physical meaning given by (1.11) only for 

Ui = U2 = Us = 0. On the other hand, the tensor F defined by (1.4) always has components given 
approximately by (1.10). In case of gravitational field, Ui = U2 = U3 = means that the source 
of gravitation is in rest with respect to the chosen coordinate system. In this system it is Fij = (jjij 
assuming that w = 0, and hence it is true in each coordinate system. 

Note that the introduced connection Sij is non-linear. The first step of approximation yields to 
the following linear connection Ffc = (r^,.) such that 5] « ^)s^si where 



Fi = 




-a2l(? 
-as/c" 
. 


1 2 





iwj, 1 c 


/ 2 




-11112/ c 




-iwz/c 

iW2 1 c 

. 


r2 = 


- 


. —iwsjc 










iwi/c 


iw2,/c 


—iwi/c 
. 


r3 = 


- 


.iw2/c 





02 /c^ 
—iwi/c 


-ai/c^ 
-a2/(? 






—iw2/c' 
iwi/c 


. 


r4 = 


- 

-iws/c 

iW2/c 


iws/c 


—iwi/c 
-a2/c^ 


—iw2/c 
iw\/c 



a2l(? 
. 



(1.14a) 



(1.146) 



(1.14c) 



(1.14d) 



and where we assumed that {Ui) = (0,0,0, 1) at the considered point. The Newton's approximation 
obtains by replacing the components F^j, for A: € {1,2,3} with zero. Using this linear connection 
the curvature tensor is easy to be calculated and in [8] it is shown that the Einstein equations of 
gravitational field 

R^3 - \gijR = SGnc-^Tij (1.15) 

are satisfied, where the energy-momentum tensor is given by T^- = pUiUj. It is also shown in [8] 
that the following equations 

(1.16) 
(1.17) 



{lJ'Fij);k + {l^Fjh)-i + {nFki)-j = 0, 

{liFij).j = AGnc-^Ji, 



analogous to the Maxwell equations, also hold at that order of approximation. Here = pUi is the 
4-current density for gravitation, analogous to the 4-current density for the electromagnetism. 



2 Some results for the geodesies according to the ILLTG 

In this section the planetary orbits in the solar system as well as the trajectories of the light rays 
will be considered. Our coordinate origin will be chosen to be at the center of the Sun. The source 
of gravitation is in rest to our coordinates system, and we can put Ui = U2 = U3 = 0. We shall 
assume that the mass of the moving body is negligible with respect to the Sun. The general equations 
of motion of a body with nonzero mass are given in section 6. We shall deal with the orthonormal 
coordinates. 
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Let fi = 1 + ^M-, where M is the mass of the Sun. Then 























GM^y 











GM ^ 






GM 






(2.1) 



and the equations of motion ^ = —Pri(l>rkPkj for i = 1, 2, 3, 4 without approximation reduce to [8] 



dt2 



GM 



(/3- 



2)x+ -^{xvx +yvy){2 



(Py _ GM 



( 2\-i/2 

where /3 = ( 1 — ^ j . Without loss of generality we assumed that Vz = ^ and that the particle 
moves in the a;i/-plane. Thus we omit the equation for i = 3. 

Remark. Using the formulas (1.1'), and the fact that /i is a scalar which does not depend from 
the coordinate system (orthonormal or ordinary), it is easy to see that in ordinary coordinates the 
equations of motion (2.2) become: 



m + 1) 
1 



(/3-2 - 2)y + -^{xv. + yvy){2 + ^ 
/3 — Inyit = Ci(= const.), 



(2.2a) 

(2.26) 
(2.2c) 



(2.2'a) 





2 d^i dX 


q d'^X 


dT-^ 


fjL dT dT 




d^Y 


2 dfx dY 


q d^V 


dT^ 


/i dT dT 




/3- 


In^ = C\(= 


const.), 



{2.2'b) 
(2.2'c) 

where ^ and ^ are given by (2.2a,b). Note that the equation (2.2c) remains unchanged. Indeed, /U 
appears to be the same in both coordinate systems and also the quotient v'^/c'^ takes the same value 
in both coordinate systems. 

Remark. Note that if we have not only a gravitational field with a potential /U, but we also have a 
3-vector of acceleration {a^, Uy, az) and a 3-vector of angular velocity {wx, Wy, Wz), then in orthonormal 
coordinates it is 



dvx 
~dt 

dVy 

~dt 

dvz 
~dt 



(2-/3 ^)ax '^{axVx+ayVy + azVz){2 + 



(2-/3 ^)ay - -^{axVx + ayVy + azVz){2 + ^gp-pyy) + "^(vzWx - VxWz) 

V 1 
(2 - (3~'^)az - -^{flxVx + ayVy + azVz){2 + ^^ ) + 2{vxWy - VyWx) 



1 


/3(/3 + 


1) 


1 






1) 


1 





) + 2{VyWz - VzWy) 



, = -— 2 K^^x + ayVy + azVz). 



(2.3a) 

(2.36) 
(2.3c) 
(2.3ci) 



The acceleration which is caused by the gravitational potential fi is included in the 3-vector of accel- 
eration {ax,ay,az). 
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Now, using that 

_d y _ Vyx - v^y 

— --arctan-- , 

from (2.2a) and (2.2b), we obtain 

d I ndip\ d^y d?x GM , , 1 

-r = —rrX -rrV = 7r^iv,,X - V-r.y)(xV^ + yv„)(2 + —— 



d^GM( 1 W/lx 

dt ^JLC^ V ^ I3{I3 + 1)) dt\r)' 



Two cases will be considered. 

(i) Assume that v « i.e. consider a planetary orbit. Then 2 + « 2, 5 and hence 

d /■ 2d(p\ _ 5 GM /■ 2d(p\ d /l^^ 
di\ 'dtJ ^ ^2~y IttJjArJ' 

by neglecting the expressions of order c~^. The solution of the previous equation is 

^ =^'^''^1^7^)' ^2 = consi. (2.4) 
Further, using the metric (rfr)^ + r^{d(fi)'^ — (P't^ = ds^ in orthonormal coordinates, we obtain 

and get 

where p = r~^. Prom (2.2c) we obtain approximately 

(l-^Y^^^ -l+c-^C 
\ (? ) r(? V 2 re? ) ^ ' 

where C3 = (?{Ci — 1) = const. Hence we obtain 

GMp G^M'^p'^ 



2 2 
V — c 



= 2Cs + 2GMp{l - - 6Cs^ - 3 



Replacing this value of v"^ into (2.5), after neglecting the summands of order c ^, we obtain the 
following differential equation 



and hence 



/, 3G^M\ dp 



where p and q are constants. Indeed, p = pi + P2, q = P1P2, where pi = P2 = and ri and 

r2 are extreme distances (ri < r2) from the planet to the center of the Sun. It is well known that 
C2 = \/ GMa{l — e^), where a = (ri +r2)/2 and e is the eccentricity of the approximate ellipse. Thus, 

/ 3GM \ dp 

\ ac [L - e )/ y'-p^ + pp - q 
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Hence we obtain 

/ 3GM \/ . p-f ^\ / 3GM \ / . 2p - pi - p2 ^ \ 

V ac^(l-e^)/v /p2 _ / v ac^(l-e^)/v pi - p2 ' 

V 4 



Assume that </5 = if p = pi, and hence C4 = 7r/2. This equation can be solved by p and its solution 
is 

Pi + P2 , Pi - pi 



H cos 



/ _ 3GM 
.'^V ~ ac2(l - £2) 



Now it follows that the angle between the points for which p = pi and p = p2 is equal to 7r+ 



6GM7r 



(2.8) 



ac^(l-e^) 



and hence the angle for the perihelion shift is equal to At^ — ^^2(1,^2) ■ 

(ii) Let us consider the trajectory of a light ray near the Sun. We denote by R the radius of the 
Sun. In this case /3 = 00, 



dt\ dtJ c2 V dt)dt\r) 



dt\ dt / c2 V dt J dt\r 

and its solution is 

2dcp ^ ( GM 
"^-^ =C2exp -2— 
dt \ rc 

Analogously to (2.5) we obtain 



r^-^ = C2exp(-2— C2 = const. (2.9) 



\dip/ 

and replacing ii = c, we get 

We shall prove now that the light has a constant velocity c in a gravitational field, using orthonor- 
mal coordinates. Summing up the equation (2.2a) multiplied by 2vx and (2.2b) multiplied by 2vy we 
obtain 

dv^ 2GM \ v^\, , , ,/„ 1 ^^;2- 

^ = ^ 1' + ^ J (^"^ + ^"^^ + (^"^ + I' + WTV)) ^- 

2GM, .V~^ 

= o-{xVa: + yVy)^ , 

1 + 0^ 

Hence, if ?;2 = at one point, then ?;2 = holds on the whole trajectory, because v"^ = satisfies 
the previous differential equation and the initial condition. 

If r = i?, then = c and from (2.9) we get C2 = iicexp(2^^). By replacing this value of C2 
into (2.10) we get 

By neglecting the expressions of order we obtain 

/dp\2 2 1 4GM 4GM 

which is equation of a hyperbola. Note that according to the Newton's theory the trajectory is a 
different hyperbola: (^)' + p2 ^ _ 2GM ^ zgm^. 
From (2.12) ip can be determined as a function of p: 

pR^^-2GM 

= arccos _ , (2.13) 
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such that (p = if p= ^. Hence from (2.13) we obtain 

2GM ( 1 2GM\ 



(2.14) 



It is easy to conclude from (2.13) that the angle of deflection of a light ray near the Sun is equal to 
This conclusion can be obtained also by the following way. Assume that the light ray is parallel 
to the X-axis at (0, i?, 0, 0). Let us put ^ = wcosa, ^ = — t;sina, x = vt, y = R and r = ^fW^^rvH^ 
in (2.2a) and (2.2b). Then 

dv 



dt 



cos a 



. da GM r 

vsma— = 



dv . 

sm a 

dt 



V cos a 



dt 

da 
dt 



GM 



V COS a , 
vt H 5 — [xvx 



/, „ vsma , 

-(l + ^)i?-^(x.. 



^yvy){2- 



1 



/?(/? + !) - 

1 ^ 

/3(/3 + l)' 



(2.15a) 



(2.156) 



Summing up the first equation multiplied by sin a and the second equation multiplied by cos a, we 
obtain 

da GM 



dt fir^ 



1 H — 2 ) (-Rcosa + vtsina). 



By putting I) = c, /U ~ 1 and R cos a + vtsma k, i?, we get 

da 2GMR 
Itt 



Aa 



-r 

J — c 



2GMRdt 4GM 



c{R'^ + cH'^yi'^ Rx? 



(2.16) 

by using orthonormal 



Although we deduced the main formulas A^p = ^c^i^li^) ^^'^ 
coordinates, they remain unchanged if we use ordinary coordinates. 

To the end of this section we shall consider the motion of a particle (with zero mass) under the 
gravitation of k gravitational objects. Indeed, we have in special case a motion in the solar system 
with perturbations of the planets and other gravitational bodies. Let rrii be the mass of the i-th 
gravitational object and be the distance from its center to the particle. Further, let the 3-vector of 
velocity of the particle is v = {vx,Vy,Vz) and the 3-vector of velocity of the i-th gravitational object 
is Ui = {uix,Uiy,Uiz) and let us put 



?7i = -cMn 1 + 



(2.17) 



Let us denote the 4- vector of velocity of the i-th object by 
iUa,Ui2, Ui3, Uii) = 



Uix ^ty Ujz Y 



such that the 4-vector potential is ^1 + {Uii,Ui2, Ui^, Ua). 



The tensors Fpq and 4>pq considered for the i-th object are 



PI 



dUj 
dx„ 



dU, 
dx, 



(do not sum over i). Now, analogously to (2.3) the equations of motion under the gravitation of the 
i-th object is given by 



dvx 
~dt 



(2-/3 '^)ax '^{axVx + ayVy + azVz){2 



/?(/?+ 1) 



) + 2{VyWz - VzWy] 



(2.18a) 
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dVy 

~dt 

dvz 
dt 



(2-/3 ^)ay ^iaa:Vx + ayVy + azVz)i2 + 



1 



m + 1) 
1 



(2-/3 )a, - -^{a^v^ + ayVy + a,v,){2 + ^ 



) + 2{VxWy - VyWx) 



d 1 



dt -iP^ c 



1 . 

--2[0'xVx + ayVy + a^Vz), 



(2.186) 
(2.18c) 
(2.18rf) 



where 



i^il^ '^'4-5 — 

0X4 OX\ 

UaT. 1^12-^ — 

0x2 ox^ 



a„ = 



I^i2-^ '^j4^ 

0X4 0X2 



Uil-^ f^iS^ 

0X3 0X1 



a, = 



Wz 



If the potential Ui is not time dependent, then 



Ox = - 



1 dUi 



? dx 



1 dU,, 



dy 



l^iS^ '^j4-5 

0X4 0X3 



Ui2-^ tVji- — 

OX\ 0x2 



1 dU, 



0^ 



^ dz 



Further, after some algebraic transformations, the final equations for the motion of the particle under 
the gravitation of all objects become 



dvx 
dt 



dUi Mvl - 2uixVx 



) + 



dx 



'-+ 



dUj MvxVy - 2uixVy dUi Mv^v^ - 2uixVz 

"To o ~r 



dy 



dz 



(2.19a) 



dVy _ \^Ej_/q-2 _ 2 ^ g "»x^x + UjyVy + UjzVz . _^ dUj Mv^Vy - 2UiyVx ^ 

dt ^ \- dy (? 



dt 



dx 



dUi MVy - 2UiyVy dUi MVyVz - 2UiyVz 

dy (? 



dz 



74- 



Er^t/j 2 _ 2 _^ r, '^ixVx+UiyVy+UizVz . ^ dUj Mv^Vz - 2UizVx ^ 
I dz 



dx 



dUi MvyVz - 2uizVy dUi Mv^ - 2uizVz 

' 7\ o ~r 



dy 



dz 



d 1 



dt Ii _ 



M/1 - ^ 



(2.196) 



(2.19c) 



1 ^ 1 fdU, dU, _^dU, \ 



where M = 2 + ^^^^jijj- It is important to emphasize that ^^Vx + ^§^'^v + ^^^^ is not a total 
differential of Ui, because Ui depends also on the coordinates of the i-th body. 
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3 Representations of the geodesies near the Schwarzschild 
radii according to the ILLTG 

In section 2 were considered the geodesies in polar coordinates with approximation, neglecting the 
expressions of order c^^ and smaller. In this section we shall consider the geodesies in polar coordinates 
without approximation, using orthonormal coordinates. Starting from the equation (section 2) 



d_ 
dt 



we obtain 



(^2dv\ _ ( od^\GM( 1 \dn\ 

V dtJ V dtJ i^c^y dt\r)' 

iWt))- 



GM / 1 \ ^ /1\ 

7^1 ^ /3(/3 + l)JdtVrJ" 

Replacing in this equation that /? = Ci + In^l + = Ci + In^l + as a function of p, we 



obtain 



s(K'"f )) = + WTT)))<?)' 



and further like in section 2 we obtain 

2 



where 



2 2 

V = c 



1- Ci+lnM 



2GM 

2 



The constant C2 ^ will be determined, assuming that = at the considered point. Indeed 



A* = 1 + 



/3(/3 + l). 
GMp 



2 2^-2 / /■ 2GAf / 1 
P^ = ^o^^2^exp(y-^(2 + ^^ 



dp 



and hence 



-T- +P =^^^exp /2 — - 2+— — 
We convenient to use the short notations 

Now, (3.1) can be written in the following forms 



1) 



dp 



GM 

■(2 + 

pc 



:2 V ^ /?(/? + 1) 



dp 



P=PO 



). (3.1) 



f^GM / 1 N , 



^(d^) + ^ _ Po 1 



\d(py 

Let us denote by / the following characteristic function 



to 




.p1 


1 ■ 


) - 




y2 





(3.2) 
(3.3) 



(3.4) 
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0,2637 0,5030 



Fig. 1 



Note that = [1 - (Ci + In ^) "^J and we can consider the function (3.4) only for those p, for which 
> 0. Two cases are possible and we shall consider them separately. 

(i) Assume that Ci > 1. The case Ci = cxd, i.e. w = c is included. Physically Ci > 1 means that if 
the particle left the gravitational filed and goes into infinity (r = oo), then it will move with velocity 

vq such that Ci = ( 1 — % ) . Specially if Ci = 1, then vq — 0. The characteristic function /(p) is 



defined for each p > and obviously limp^o f{p) — 0. Since 2 < 2 + -g^^q^ < 2^, we get /i"* < < /i^. 
Thus, if p ^ oo, then v ^ c and using the previous inequalities we obtain limp^oo /(p) = 0. We shall 
find the maximal value of /(p). Indeed, f'{p) = can easily be written in an equivalent form 



dv vGM / 1 
Tp'^ ixc^ \ /?(/?+ 1) 



p2 = vp. 



By differentiating the equality /3 - In p = Ci, it is easy to get ^ = (^1 - Thus, /'(p) = 

if and only if 

GMpl/ w2\3/2 GA/p/ 1 
~i^vV~~^) ^'"^I^V^ [3{f3 + l), 
and after some algebraic transformations it it is equivalent to 



,2 



GMp. (3.5) 



The physical meaning of (3.5) is very simple. Indeed, it is easy to verify that the equations (2.2) 
are verified if a particle moves with velocity v — y/GMp on a distance r = ^ on a circle. It is just like 
in the Newton's theory, but remember that Ci > 1 which is analogous with ^v"^ — GMp > in the 
Newton's theory and according to it the trajectory is hyperbola, but not a circle. We shall prove that 
this effect in ILLTG is a " relativistic effect". The equation (3.5) is equivalent to 



1 . GMp^ 



and will be called Schwarzschild 



The graphic of the right side of (3.6) is given on Fig. 1. 

Each of these circles has radius in the interval -^jj , ^ ^^^^ ^^^^ 

radius which corresponds to the chosen value Ci > 1. We shall denote these radii by index s. The 
graphic of the function /(p) is given on figure 2. On fig. 3 is given the graphic of /(p) when Ci = 1. 

2 2 

The graphics are compared with the corresponding Newton's case /Af(p) = ^ = const +2GMp • 
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Fig.2 



P 




Now according to (3.2) we are able to draw some conclusions for the trajectory. Let be given p, 

V and ^ at a chosen point. Since p and v are known, the constant Ci (> 1) is known and also the 
function f(p) is known. We have the following options and the conclusions follow by Propositions 1 
and 2 at the end of this section. 

^ dp \2 I ^2 

- If "^^^2^1 — > f{Ps), then the equation (3.2) has no solutions for po a-nd hence ^ = will never 
be achieved and the particle will tend to the center r = 0. This easily can happen for large initial 
values of 

dip 

- If "^^^2^1 — — f{Ps), then the particle will move on a circle. 

- If ''^2eJ — < fiPs), and p < ps, i.e. r > r^, then the trajectory of the particle will be out of the 
Schwarzschild sphere (S.s.), it will look like a hyperbola and = will be achieved once. 

- If ''^^^j — < f{ps), but p > Ps, i.e. r < Vg, then the trajectory of the particle will be only inside 
of the S.S., and it will tend to the center r — 0. The condition ^ — Q will be achieved once. 

According to the previous discussion we see that the S.s. with radius r = is a sphere of 
separation. We shall consider further in more details the motion of a light ray. In this case Ci = oo, 

2 

V = c and ps = -^jj- On fig- 4 in the xy-plane are given different trajectories of a light ray assuming 
that ^ = at the initial point. 

If the light ray cut the S.s., then it will tend to the center r = 0, because the equation (3.2) has 
no solution. The equations of motion of a light ray are simpler, because /3 — oo. The equation (3.3) 
for light rays becomes 



/dp\2 ^ r pj p' 1 / ^ GMp Y 



V dip 



-) , (3.7) 
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Fig.4 



while the equation 

dt 

becomes 

2^'^ r< (a ^ GMpy 



Assuming that ^ = if r = tq, wc find 
and hence 

, I I GM „ 

The total angle of deflection of a light ray outside the S.s. is equal to 

Aa = 2 , , "'^ = - TV. (3.9) 

(ii) Assume that Ci < 1. PhysicaUy it means that the particle can not leave the gravitational field 
because of small energy. In this case f{p) is defined only for p > p^,, where is the solution of the 
equation 

Note that limp_+p»+ f{p) =00 and limp^oo f{p) — 0. Similarly to the case (i), the condition |^ = 
is equivalent to (3.5) and further it reduces to (3.6), but now Ci < 1. Now three subcases are possible. 
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(iia) Assume that 0,9313 < Ci < 1. Then the equation (3.6) has two solutions (fig. 1), pc and 
Ps- The smaller solution pc has just the same meaning as in the Newton's theory: Each particle on 
a distance Tc = — will move on a circle with velocity JGMpc.. The larger solution p^ has the same 

Pc 

property but it is a relativistic effect, because ps > 0,2637^$^-, (r^ < 3, 7921^p^), and fs — y will 
be called again Schwarzschild radius because it separates the curves just like on fig. 4. The graphic 
of the function f{p) is presented on fig. 5. 

Analogously as in the previous cases we have the following conclusions. 

^ dp ^2 I ^2 

- If '''^^gj — > f{Ps), then the equation (3.2) has no solutions, such that on the trajectory will 
never be ^ = 0. Thus, the particle will tend to the center r — 0. This easily can be achieved if we 
throw a particle with small velocity which is almost radial motion, and then the particle will come to 
the center of the Earth, assuming that mass of the Earth is concentrated at a point. According to the 
Newton's theory the particle will move on an ellipse with large eccentricity. 

^ dp j2 I ^2 

- If ''^2eJ — — fiPs), then two cases are possible. If p — p^, then the particle will move on a circle 
on the S.S., while if p < ps, then the particle will tend to the center r ~ 0. 

^ dp n2 I ^2 

- If f{pc) < '''^^gj — < fiPs), then two cases are possible. If p < ps, then the equation (3.3) 
has two solutions smaller than ps and then the trajectory will be like an ellipse, analogously to the 
planetary orbits. The distance r to the center of the gravitation satisfies ri < r < r2, where pi — ^ 
and P2 = ^ can be found as solutions of (3.3). On the other hand, if p > ps, then the particle will 
remain in the S.s. and will tend to the center r — Q. 

- If '''^agj — — f{Pc), then two cases are possible. If p < ps, then the particle will move on a circle, 
while if p > Ps, the particle will remain in the S.s. and will tend to the center r = 0. 

- If ''^y2^i — < f{Pc), then the particle must remain in the S.s. and will tend to the center r = 0. 
{lib) Assume that Ci = 0, 9313. Then the equation (3.6) has a double solution p^ = pc = 

0, 2637jTp-. In this case ps and pc coincide, the graphic of /(p) is given on fig. 6. The discussion is 
analogous as previous. 

(iic) Assume that Ci < 0, 9313. Then the equation (3.2) has no solution and the particle tends to 
the center r = 0. The graphic of /(p) is given on fig. 7. 

At the end of this section we give some global conclusions in condensed form. First note that the 
Schwarzschild radius is defined for each constant Ci > 0, 9313, belongs to the following segment 

GM 1 CM GM 

<rs< TT- = 3, 7921^5- 

c2 - ' - 0, 2637 c2 ' c2 

and Ts for the light rays is equal to 

- If at the initial moment the particle is inside the S.s. and the equation (3.2) has a solution p > p^, 
then the particle has never been outside the S.s., it will remain in the S.s. and will tend to the center 
r = 0. 
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y 




- Suppose that at the initial moment the particle is outside the S.s. If there exists an unique 
solution p of (3.2), such that p < ps, then the trajectory will look like a parabola or hyperbola and 
will never pass inside the S.s. If there exist two solutions pi, p2, {pi,P2 > Ps), then the motion will 
look like an ellipse. 

- If the equation (3.2) has no solution, then it will never be ^ = on the trajectory. The particle 
will tend to the center r = 0. In this case the particle can cut the S.s. 

We shall explain why the other cases are not possible. If the equation (3.2) has at least one 
solution, the trajectory will be axially symmetric with respect to the line which connects the center 
r = and the point where = 0. Moreover, the trajectories are "time reversible", which means that 
if a particle moves in one direction on the trajectory, then it can move also in the opposite direction 
of that trajectory. 

We shall prove the following statements. 

Proposition 1. // a particle cut the S.s., then inside the S.s. can never be ^ — 0, and the 
particle tends to the center of gravitation. 

Proof. Assume that p, v and ^ are given. Then p and v determine the value of Ci. Now the 
equation (3.2) can not have a solution, because if such a solution p' exists, then p' < p, because f{p) 
is a decreasing function for p > ps. This finishes the proof. 

Corollary 1. The trajectories of "elliptic type" can not cut the S.s. 

Proposition 2. Inside the S.s. there does not exist a trajectory of "elliptic type". 

Proof. Assume that there is a trajectory of elliptic type which is inside the S.s. and let the extreme 
distances of the particle to the center of gravitation be ri = ^ and r2 = ^ . Then there will exist a 
radius r G (?'i, ^2), such that for the same value of Ci, a new particle can move on a circle of radius 
r. This is a contradiction with the fact that each S.s. is not contained inside another S.s. because the 
radius of the S.s. is minimal such that a particle can move on a circle for the same value of Ci. 
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4 Some relationships between the geodesies in ILLTG and GR 



In this section we shall find some connections between the planetary orbits according to ILLTG using 
orthonormal coordinates and GR using the Schwarzschild solution. Namely, we shall show that there 
are simple transformations which the geodesies from ILLTG map into the geodesies from GR and 
conversely. 

Let us denote hy r, p = ^, (p, t the polar coordinates in ILLTG, and by r*, p* = 93,, t» we 
shall denote the polar coordinates in GR. Note that in both tests in the solar system (deflection of 
the light ray near the Sun and the perihelion shift) were measured angles and these angles were the 
same for ILLTG and GR. This associates us to assume that (/j, = ip, i.e. we can choose the angle ip 
as a common parameter, and then we can find the relationships between p and p*, as well as between 
t and t*. These relationships will be found by neglecting the expressions of order c~^. We convenient 
to write (p instead of . Although we shall consider ellipses as planetary trajectories, the obtained 
results are also true if the trajectory of a planet is a hyperbola or a parabola. 

Let us start with the equations 



d(p = 



dip = 

dp.. 



(1 



> GM 

'ac2(l-e2) 



)dp 



Vipi - p*)ip* - P2) 



■ p){p- P2) 

^ GM 

1 H T^P* 



2GM 



ac2(l -e2) 



from the ILLTG and GR respectively. The solution of (4.1) we found to be 



Pi +P2 _^ Pi_ 



P2 



■ cos 



3GM 



/ 6UM \ 



ac2(l-e2). 



(4.1) 
(4.2) 

(4.3) 



and the solution of (4.2) is given by 



-(> 



2GM 



)( 



+- 



GM 



■(- 



ac2(l-e2)y V 

- P*){P* 
3GM 



arcsm ■ 



2/0* - Pi- P2 



Pi 



P2) + 



Pl + P2 



ac2(l - e2) 



)( 



arcsm ■ 



. (- 

2p* - Pl- P2 



arcsm - 



P2 

2p* 



Pl - P2 



Pl - P2 



2) 



Pl 

GM 



P2 



+ 



I) 



\/{pi- p*){p* - P2) 



and hence 



cos 



^(1 



3GM 



ac2(l -e2). 



\ 5-V (Pi 

pl - P2 & ^ 



P*){P* - P2)\ 



'2p^ - pl - 



pl - P2 



I.e. 



cos 



(- 



3GM 



2p* - Pl 



ac2(l-e2)yj P1-P2 
By eliminating the angle (p from (4.3) and (4.4), we get 

GM 



P2 GM2(pi -p,)(p. 



P2) 



Pl - P2 



P* + ^-(Pl - P*)(P* - P2) 



and also the converse relation 



p* 



GM 



{pl- P){P- P2)- 



(4.4) 

(4.5) 
(4.50 
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Now we shall find the analogous equation of (4.1) written in ordinary coordinates. We shall use 
the letter R for radial distance in ordinary coordinates. Fist let us write the solution (4.3) in the 
following form 



Pi+ P2 , Pi - P2 



■ COS 



3GM 



ac2(l-e2) 



and after differentiating, using that dr = ^1 — ^^^dR w ^1 — ^^^dR, we obtain 

fi 



Pl+ P2 , Pi - P2 
1 COS 



After integration of this equation we get 

'Pi + P2 , Pi - P2 



GM R 

^ 2- -W 

&■ Ri 



COS 



^ 1 



3GM 



ac2(l - e2) 



2 2 

where the constant of integration is chosen such that R = Ri = (pi)"^ for = 0. Further, for elliptic 

Ri ~ Ri 

3GM Nm-l 



trajectories with small eccentricity we can put In « — 1, and hence 



Pl+ P2 , Pi - P2 

— 1 — cos 



(4.6) 



GM 



ac2(l-e2) 

This function for = and </9 = f (^+ aJ(i-e)i) ) takes values (pi)"^ and (p2)"^ - 
(P2)"^(l + ^M£p£i)^ ^ respectively. Thus, in (4.6) p2 should be replaced by ^2(1 - ^MfpSil^ 



Hence we obtain 



PI+P2-P2 



GM{p2-pi) 



P1-P2+P2 



GM{p2 - pi) 



cos ip ( 1 



3GM 



ac2(l - £2) 



After some calculations, one can verify that the function R ^ satisfies the equality 

GM 



R-' 



-{pi- P){P- P2), 



I.e. 



p. = R-\ 

The differential equation for R~^, analogous to (4.1) is given by 



dip 



4 



V(pi - R-'){R-' - P2) 
Note that if we put r = r^[l- ^) = r 

\/{Pl* - P*){P* - P2*) 



, GM 1 

1+ ^^T^ + 



2GM 



d(fi = — 



c2 R ac2(l-e2)J 

then the equation (4.1) becomes 
2GM 



(4.7) 
(4.8) 

(4.9) 



GM 
1 H j-p* 



ac2(l-e2) 



(4.90 



where pi* = (n + and p2* = (ra + ^)-^ 

The equation (4.8) shows that for v << c, the trajectory of motion in the (r, tp) plane is 
the same with the trajectory according ILLTG using ordinary coordinates. 

Further we shall find the relationships between t, and T. According to GR it is 



* dt* V 



1 



2GM 



-p* 
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^dif ^ 2GM 1 GMp\ 



and according to ILLTG it is 

= (^2(1- 

The quotient of these two equations yields to 

\r) dU H "*" /?(/? + !) c2 
Since ^ = 1 + ^(pi - p){p- P2) according to (4.5'), and J(Jlp^ » 5 we obtain 

dt ^ 1 GMp GM , ,\ 

Moreover, C3 = 1, because in case of parabola or hyperbola, ^ > 1 if p ^ 0. Hence we get 



dt , IGMp MM, 



Assuming that the trajectory is elliptic with eccentricity close to 0, then 

dt _ 1 GMp 

Moreover, we can find the connection between t* and the ordinary time coordinate T from ILLTG 
using that dt = [l ^ ^§^^dT, i.e. 

dT _^ 1 GMp 
dt» 2 c 



(4.10) 



The formula (4.10) docs not show how much the time according to GR is "faster" (or "slower") 
than the time according to ILLTG, but only a factor for transition from the geodesic of GR into 
geodesic of ILLTG in ordinary coordinates. 



Further we shall compare also the null geodesies, i.e. the trajectories of the light rays. According 
to (2.14) in ILLTG it is 

1 GM 



p = —cosip + -f^^{2 - 2costp), (4-11) 



while according to GR it is [5] 



Hence we obtain 



1 GM , , , , , 

p* = -cos(^+— ^(2-cos>)- (4.12) 



I.e. 







GM 


p* 


= P + 


- ) 

i?2c2 


p* 


= 1 + 


GM 


p 


pR^c^ 


r* 




GMr 


r 


= 1 - 


R^c^ 



cos <f{2 — cos (f), (4-13) 
cos<^(2 — cos<p), 



cos 1/3(2 — cosy). (4-14) 



Further, using (4.14) as well as the following two equalities 
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which are almost the same in both theories, we get 

2 dt 



fr*Y dt 
\ r J dt^ 



dt ^ GMr ^ , 

^ = 1 + 2^^ cos (^(2- cos (p). (4.15) 

The angle of deflection from a point A to a point B when the ray passes near the Sun according 
to (2.16) in ILLTG is equal to 




= ^ 2W1 - _ +2W1 - , (4.16) 



and also using PPN coordinates GR yields to the same formula (4.16). This value (4.16) is confirmed 
with a very precise measurements in the recent test with Cassini spacecraft [4]. According to ILLTG, 
the angle of deflection is always twice larger than the classical Newton's result, and can be explained 
with the coefficient 2 — = 1 + /(? = 2 in (2.2) when v = c. Note that this deflection of the 
light is not a "gravitational effect", because it will appear if we have only an acceleration without 
gravitation, and ILLTG differs slightly in these two cases when there is a gravitation and there is only 
simulation of gravitation using inertial forces. About this will be said more at the end of this section. 

We shall compare now the radial motions of a particle with large velocity ii in a gravitational field. 
The Newton's relation — GMp = const, according to ILLTG should be replaced by (2.2c), i.e. 

GMp 



- ln(l + = const. (4.17) 



According to GR, for a radial motion in a gravitational field we get 



_ 2GMp 



7T- 



const. (4.18) 



where c, = c(l — ^^^-^) is the velocity of the light measured in the gravitational field. The physical 
interpretation for both formulas (4.17) and (4.18) obtains if we multiply them by moc^. Prom (4.17) 
we obtain 

Ek+Ep = const. (4.19) 

because in section 6 will be proved that the potential energy for a particle with small mass is equal 
to Ep = — moc^ In^l + • The equation (4.18) can be written as 



Etotai y 1 - = const. (4.20) 

If is very small, i.e. we consider a weak gravitational field, (4.18) can be written in the form 

1 ^GMp ^ 

C — 5 h 0, C = const. 



and hence 

d 1 _ GMdl 



dt fi ^ dtr 
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I.e. 

a 



, = C^v, (4.21) 

where a = is the acceleration and v = dr/dt. The analogous equation from (4.17) in weak 
gravitational field is the following 

d 1 a , , 

= = 3^- (4-22) 



dt 



In a short time interval the acceleration a is almost constant and the equations (4.21) and (4.22) 
should be close to the equation of motion of a particle in STR under a constant acceleration a. The 
STR motion 



is just the equation (4.22). Indeed, in (4.23) a is a constant, because the acceleration a is unchanged 
under the Lorentz transformations of motion in the x direction, because the tensor with components 



- 








a' 


























.—a 








0. 



is unchanged under such Lorentz transformations, analogously as the component of the electric 
field is constant under such Lorentz transformations. 

Both of equations (4.21) and (4.22) in the (a;t)-plane are hyperbolas, but with different param- 
eters. The parameter C has no physical moaning according to STR for a motion under a constant 
acceleration, although this constant in (4.18) has a good physical meaning. Thus, the equation of 
radial motion is not close to the STR. Moreover, from (4.18) one can conclude that ^ a(l - 3^), 
such that the acceleration changes its sign if v > c/ \/3. 

Wc note that ILLTG is very close to the STR. Namely, if wc consider a nonincrtial system without 
gravitation but there is only a 3- vector of acceleration {ax,ay,az) and a 3- vector of angular velocity 
{wxjWyjWz), then in the equations of motion (2.3) we should only put /i = 1. This is obtained in 
the frame of STR, under the assumption that the Jacobi matrices with respect to the incrtial systems 
are orthogonal and hence the coordinates in such a system are unholonomic [9] . The influence of the 
coefficient fi in (2.2) and (2.3) is very small. For example, if we have a noninertial system without 

gravitation with the same acceleration {ax,ay,az) = (^—^^x,—^^y,—^-z^, then the angle of 
deflection of a light ray near the Sun will be again Aa = , while the angle between two perihelions 
will be Aa acHi%) ^ instead of Aa = ^^3^. The equation (2.4) wiU not be changed, and the 
equation (4.17) will be replaced by 

1 GMp 

= — = const. 



5 Two proposal tests 

After the discussion in section 4 we propose the following tests for which we have different results 
from ILLTG and OR. 

Test 1. Let Ai, A2 and ^3 are three not close points of a planet orbit or a spacecraft trajectory, 
and let LAqOAx = (po, LA^OA-i = t^i, LAqOAj, = (^2, where Aq is the perihelion of the planet orbit 
and O is the center of the Sun. We denote by Ti and T2 the times for which the planet comes from 
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Ai to A2 and from A2 to A3 respectively. We will calculate Ti : T2 according to ILLTG and GR, and 

later can be compared with the experimental results. We shall find Ti : T2 where Ti and T2 are the 
proper times for the Earth, where the trajectory of the Earth is observed. Let us denote the proper 
time for the Earth by r. 

First let us calculate Ti : T2 according to ILLTG. Prom (2.4) we obtain 



2# odif dt dT ^/ 5GM\/ GM\/ GM\ ^/ 1GM\ 

dip 1 



Thus 



and hence 

where p is given by (4.3). Using that is close to 0, and replacing t^^l — ^^J^j^sy^ with a new 
variable, one obtains 

T T = ^ C(/9-^(v?i)v?i - p~'(yo)yo) - A{p-^{^2)^2 - p-^(yi)yi) 

C ac2(l-e2) G2 ^ 



where 



pifi pLfi pip2 r^2 

A= p~'^dip, B = p'^dtp, C = p~^dip, D = I p'^dtp. 
After solving these integrals, where p= \ [{pi + P2) + {pi — P2) cos<p], we obtain 

_ Pl+ P2 ^/prp2(cot ^ - cot ^) P2- Pl ( cot ^f- cot ^ 



■ arctan ^ ■ + 



( ^"^2 2 A (r^. 

\0■> + 0^ COt^ ^ + rti cot^ ^ ' ^ 



(P1P2)^/^ P2 + Pi cot ^ cot ^ P1P2 Vp2 + Pi COt^ ^ P2 + Pi COt^ ^ - 

Pi + P2 Vplp^(c0t ^ - cot ^) P2 - Pl / cot ^ COt^ N 

(piP2)3/2 P2 + Pl cot f cot + P1P2 lp2 + Pl COt^ ^ p2 + pi COt^ ^ ' ^''•"'^ 



2 -v/PiP2(tan ^ — tan ^) 

J= arctan ^ , , 

y'piP2 Pl + P2 tan ^ tan ^ 

2 ^pip2(tan ^ - tan ^ 

^= arctan ■ — 

VP1P2 Pl + P2 tan ^ tan ^ 

Now let us calculate Ti : T2 according to GR. Using that = Gi, and (4.5') we obtain 

d(p dip 



D=-^ arctan Vp^p^v- ; ■ (5-5) 



P^Gl (^^^2 



1 



^(^-1)(1-^ 
p p 



dip, 



where p is given by (4.3). If the eccentricity of the ellipse is small, then the expression ^^^^(y — 
1)(1 — y) can be neglected. Then, analogously to (5.1) we obtain 

A 3GM G{p-\ipi)ipi - p-^ipo)ipo) - A{p-^ip2)ip2 - P"'(yi)yi) 
^''■^'=C- acHl-e^) G^ ^^■^> 
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where A, C, B and D are given by (5.2-5). 

It is convenient to assume that ipo = 0, (pi = n/2 and ip2 = tt. Then 




A= arctanj^+ , '^^ , C= ^ '.^t arctanj^- , f , (5.7) 

" (P1+P2)P1P2 

, , , (5.8) 

V Pi \/PlP2 

Remember that both formulas (5.1) and (5.6) are obtained for the observers on the planet whose 
trajectory is observed. Indeed, it means that these formulas are valid only for the trajectory of the 

Earth, because wc observe from the Earth. Further let us modify these two formulas assuming that 
the observer is far from gravitation in the solar system. Then, according to ILLTG we obtain 

d(^(l + 2^) 1 hGMp- 



rfy,(l + 2^) 1 _2ri^5GAfpx 

T,:T.^ (p- + -^-^y^ ■■ I (p-^ + -2^y^' 



and hence 



^ C ac2(l-e2) C2 ^ 

f,GMBC-AD 

where A, B, G and D were previously determined. Analogously, according to GR instead of (5.6) we 
have the formula 

T T - ^ C{p-^{^i)^i - p-'(y^o)y^o) - A{p-^{>fi2)>fi2 - p'^j^iWi) 

' ■ ^ C ac2(l - e2) C2 

GMBC-AD 

+2^2 (5-6') 

where A, G , B and D were previously determined. 

Finally, if we measure from the Earth the trajectory of another planet, then we should modify the 
formulas (5.1') and (5.6') because the trajectory of the Earth is not a circle, the "speed" of the time 
is proportional to (1 — ^^)^, where r is the distance from the Earth to the center of the Sun. Thus, 
let us denote by vi and V2 the average values of the " speed" of the time on the Earth during the first 
and the second time interval. Thus, on the Earth we shall measure the quotient 

^ = (5.9) 
U2T2 y2T2 ^ ' 

where Ti : T2 is given by (5.1') and (5.6') respectively according to ILLTG and GR. 

Note that in this test we should measure only the quotient of long time intervals, large angles and 
we need also the quotient pi : /92 = (1 + e) : (1 — e). The perturbations from the other planets. Moon, 
stars, galaxy and universe (see sec. 8 and Ref. [10,11]) are not included in the final results. 

Test 2. Let us consider the second Kepler's law, i.e. consider the expression on the trajectory, 
where t is the proper time. We shall present the formulas assuming that we observe the trajectory of 
the Earth from the Earth. 

According to ILLTG we have 

^ d^d^dT ^ / 5 GMx / , GMw, , GM^ / _ 1_GM^ 
dT dt dT dr V 2 rc^ J\ rc-^ J\ rc^ J V 2 rc^ J 



22 



and according to GR (omitting the star indices), it is 



r 



C", 



where C and C" are constants close to C = ■^/GMa{l — e^). Let P and Q be two points on the orbit 
and let us denote the corresponding distances to the center of the Sun by ri and r2 respectively. For 
example, P can be the perihelion of the Earth and Q aphelion. Using the previous formulas according 
to ILLTG we obtain 



„2 



Igjl l gMfa-n) ^ eMa(l-,^) (5.10) 

\dT/Q \dT/p V 2 rac^ 2ric"'/ 2 nrac^ 

0. (5.11) 



ydrJQ \d.T ) p V 2 r2C^ 2 nc^ J 2 rirac 
and according to GR we obtain 



'\dT J Q ' ^ydrJ p 

The perturbations from the other gravitational objects are not included in (5.10) and (5.11). 
Note that in the tests 1 and 2 the departure between the GR and ILLTG is proportional to i ^ 

11 

The coefficient ^ here is just ^^^^pjy- 



6 Mass and energy in gravitational field according to ILLTG 

If a particle with mass m, measured when it is far from gravitational field, has come in a gravitational 
field with potential /i = 1 + ^j^, then its mass measured from an observer far from gravitation, 
approximately is equal to j^jtm-- We convenient in this section that M,mi,m2, ... will mean proper 

mass, measured when both the observer and the mass are far from gravitation. 

Let us consider two particles with masses mi and m2. Assume that the distance between their 
centers is R. Then the gravitational force between them is given approximately by 

mi mz 1 , , 

^-"^l + ^l + ^i?- (6.1) 

because each of the masses changes under the gravitational field of the other particle. This formula 
becomes more precise if one of the masses is very close to zero, but (6.1) is of sufficient accuracy. For 

GmM 

example, if mi « is mass of a small particle and m2 = M is the mass of the Sun, then / = J^aM 

GM 

and the acceleration for the particle is a = /^J^ , which is in accordance with the formulas in section 
2. 

The energy which needs to distinguish both objects on a large distance is equal to 



E = I fdr = I G r, -dr. 

/ij 1 + i 1 + 2 



mi — m2 

which corresponds to the defect of mass 

mim2 
Am = 



/ fdr = 
Jr Jr 

, / Gmi\ , / Gm2M ,„ 



E = 



mi — m2 



If wc assume in (6.2) that m,i = M is a mass of a largo body for example a star, and m,2 = m, has 
negligible mass compared with M, then according to (6.3) for the potential energy of the particle we 
get 

E, = -mcHn(l + ^). (6.4) 
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Later we will return to this formula. In opposite case, if mi = m2 = m, analogously to (6.2) the 
potential energy is equal to 

= (6.5) 
and the corresponding defect of mass is equal to 



This energy can be written in the form 



mc 



2 



^p = T-g;^-^C' (6-7) 

^ Rc^ 

which is analogous to the formula for the kinetic energy according to the STR. 

Further, let us consider the case when i? 0. Namely, we shall consider an "ideal fusion", 
assuming that both particles are concentrated into points and assuming that R tends to 0. From (6.1) 
it just follows that the force tends to 



f=-Q, (6.8) 
while the acceleration for each particle is of type 

a ~ ^. (6.9) 

Note that this acceleration docs not depend on each of the masses. Probably it docs not dc;pcnd 
on the nature of these particles, for example their electricity if they are charged. In this case, when 
i? — > 0, according to (6.2) for the defect of the mass we obtain 

Am=^^^ln!^, (6.10) 
mi — m2 m2 



such that the mass after an ideal fusion is equal to 



mim2 mi 

m = mi + m2 m — . (6.11) 

mi — m2 m2 

Using the differential calculus, it is easy to obtain the following inequalities: 

mim2 . mi mi+m2 , , 

m2 < m — < < mi, (6.12) 

mi — m2 m2 2 

assuming that m2 < mi. Specially, if mi = m2 = m, then according to (6.6), the defect of mass is 
equal to Am = m. Hence the conclusion that in "ideal fusion" we get the largest quantity of energy 

if mi = m,2 = m. Note also that if m2 ^ in (6.12), then Jmilis- in I2± ^ 0. 

^ ^ \ n mi— 7712 7712 

Analogously to the inequalities (6.12), in general case one can prove the following inequalities 

m2 ^ mim2 1 + 



1 + Gmi - TOi - m2 1 + 



1 

< - 

- 2 



f ^1 \ f 
, ^1 r' 

^ 1 _(_ '^7712 



if m2 < mi. Indeed, in special case if we put r = 0, we obtain the inequalities (6.12). 
According to (6.4) the equation (4.17) can be written as 

Ek + Ep = const. (6.14) 
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where Ek is the kinetic energy. The equation (6.4) was obtained under the assumption that the mass 

of the particle is almost 0. Now. if wc want to obtain the equations of motion using that the mass of 
the particle is not almost 0, for example motion of binary stars, wc should use the potential 



which is obtained from (6.2), where m and M are the masses of the two bodies. The "correcting 
multiplier" 1 + appears if we want to use this potential to derive the equations of motion for the 
body with mass m. Indeed, the body with nonzero mass m influences on the gravitational potential 
of the other body. This potential (6.15) is sufficient to derive the approximation of the first post- 
Newtonian order (sec. 7). At the end of section 7 will be proved that the "correcting multiplier" must 
be 1 + 

If we consider the motion of a body with nonzero mass, under the gravitation of k bodies, then 
we can use the equations (2.19a,b,c,d), but the potential Ui instead of (2.17) is given by 



GM 



Ui 



= - (1 + 7^ j 1- (6-16) 



Gnii 



where m is the mass of the moving test body. 

Note that (2.19) does not consider the interaction among the given k bodies, but only assumes 
that their positions and velocities are given. But it docs not matter, because if it is given a system 
of n bodies with masses mi, TO2, m„, with initial positions and initial velocities, then we can 
find the coordinates of all bodies as functions of t by considering a system of 3n + n equations, i.e. 4 
equations (2.19a,b,C;d) for each body, where the potential for each pair of bodies is given analogously 
to (6.16). Indeed, then we have the following system 

(Pxj _ [ dUjj / 2 g 2 / dxjdxj dyjdyj dz, 1 dUjj / dxj 2 dxjdxj 

dt^ ~ . ^.\dxi c'^K dt dt dt dt dt dt J J dxj \ ^ dt ' dt dt 



1 dUi^ ( dxi dy^ dxi dyj\ 1 dUa 1 dxi dzi dxi dzj\i ^ 



dy_ 

dt^ ~ A dyj c^Kdtdtdtdtdtdtj) ci^ dxj \ ' dt dt dt dt 



-A y dUjj / 2 Q, 2 / dxjdxj dyjdyj dzidzj \\ 1 dUjj / dxj dzj dzidxj \ 
dt^ .^.\ dzj c2 Wt dt dt dt dt dt )) dxj V dt dt dt dt ) 



where {xj,yj,Zj) are the coordinates of the j-th body, /3j = ^1 — ^(j,-^)^ + i'si')^ + ("^)' 



-1/2 



and 

c'^rrii , -'■ + 7~^T-r/, 3Gto 
rrii -rrij 1 + iilli^ ^-^ \ rsjC 



'^--;;r^-^^^n(l + 7^)' r-,, = ^(a., - x,)^ + (y, - y,)^ + (., - .,)^. (6.18) 
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Note that + + is not a total differential like in (2.19d). 

In the general equations of motion (6.17) are not included the multiple moments, the gravitational 
radiation and the decrease of the orbit periods of the binary systems. Indeed, these equations lead 
only to the first post-Newtonian order, because of the weak choice of the potentials (6.18). But, using 
the ordinary coordinates we have the same eakudations like in GR about the gravitational radiation 
and the decrease of the energy and orbit period of the binary systems. 

At the end of this section we shall make some considerations about the atoms of hydrogen and 
deuterium. First we accept the following convention. If we have a particle with mass m and electricity 
e, then we convert this electricity into mass ±-y=== (the sign depends on the charge positive or 

negative), such that the Coulon's force / = j^^^ can be written as a gravitational force / = 
Q mim2 _ -pj^jg ijia,gg ^iii be called "charged mass". In order to distinguish the charged mass from the 
ordinary mass we can write m ± z -^^7== , where i is the imaginary unit. According to this identification, 
the mass of the proton is equal to 

rup = (1, 672415 x 10"^'' + i x 1, 8594077 x 10-'^)g. 

Let us consider the fusion of two atoms of hydrogen into atom of deuterium. According to the 
previous discussion, we will obtain the energy 

where e = 1, 8594077 x 10~^g. If we assume that R = 2x 10~^'^cm is the distance between the centers 
of the hydrogen cores just before the fusion, we get energy E = 2,5670184 x 10~^^ gxc^. On the 
other hand, it is known from the nuclear physics that this energy is about E = 2, 569 x 10"^'' gxc^. 
Further let us consider the question for calculation of the coefficient n inside the atom. First 

we should remember that in case of gravitation it is = 1 + and remember that = v'^ if 

2 

the particle moves on a circle with velocity v. Hence it is better to write /i = 1 + Now, let us 
consider the potential fi, when an electron is moving around the core of the atom. If its velocity is 
approximately v = c/100, then /x « 1 + 100^^ = 1,0001. If we consider another particle instead of 
electron, for example a non charged particle, since the gravitational force is much smaller than the 
Coulon's force, its velocity around the core will be much smaller and hence ji will be almost 1. Hence 
the conclusion that the potential /i depends on the nature of the particle, i.e. it is not the same for all 
particles as it was in the previous sections. 

Let us calculate the potential /i for the proton near the core of the hydrogen on a distance 2 x 
lO^^'^cm. We can not use the formula = 1 + where e is charged mass, but should use the 
following revised formula 

M = 1 + ;^A, (6.19) 

where the coefficient A should be such that = v'^. It is easy to see that the coefficient A should be 
the quotient between the charged mass and the gravitational mass, i.e. 

1,8594077x 10-^^ 
1,672415 x 10-24 

Hence, if we replace in (6.19) that e — 1,8594077 x 10-^g, because the ordinary mass is negligible, 
and i? = 2 X IQ-^^cm, we get M = 1, 0007674584. 

If we calculate the coefficient n for the electron close to the core, then A will be 1836,12 times 
larger than the previous, because the mass of the electron is 1836,12 times smaller than the mass of the 
proton, while they have the same charge. In this case, for an electron on distance 2, 81794 x 10-^'^cm, 
i.e. radius of the electron, we obtain /z = 2 and /j. = 2 achieves just for the Schwarzschild radius for 
the light ray (sec. 3). 
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Now, let us return to the potential /U for a proton near the core of the hydrogen. Since the both 
masses are equal, the energy obtained by fusion is equal to 

2fm - — = 2 X 1,672415 x lO'^^g x ^c^ = 2,56505 x lO'^^g x c^. 

7 Periastron shift of the binary systems 

In this section will be calculated the periastron shift of arbitrary binary system. For example, we shall 
assume that it is a binary system of a pulsar and its companion. Applying the equations of motion 
(6.17), or simply the equations (2.18) or (2.19) with the potential (6.15) or (6.18), we shall obtain the 
required periastron shift. 

Lot m be the mass of the pulsar and M bo the mass of its companion. Let us choose a coordi- 
nate system having origin at their barycenter, in the moment when they are at the periastron. Let 
{x, y, 0) be the coordinates of the pulsar, and let {x*,y*,0) be the coordinates of its companion. The 
gravitational potential at {x, y, 0) caused by the companion is given by (sec. 6) 



U 



Rc' ) M-m 1 + ^ ' 



GM 



where R = \/{x — x*Y + [y — y*Y- Since the following potentials 

\ Rc^ ) M + m \ Rc^ J ^ ^ 

and 

GM l G^M{M-5m) 
R ^ 2 R-'c' ' ^ ' 

do not differ from U up to c~^, the same result obtains up to if we use alternatively the potential 

U' or U". 

Let us denote the 4- vector of velocity of the pulsar by 



(y.) = ^^(^^,^,i), 

/-|^ _ 1,2 V IC IC ic ) 



and the 4- vector of velocity of its companion by 

(f/0 = ^^fe^,^,i). 

/i u2 V IC IC IC ) 

Using the potential (7.1') we can calculate the components ax,ay,az,Wx,Wy,Wz. which appear in 
(2.18). It is convenient to use the notation r = \J x^ -Vy"^ and r « M+m ^' cosa = ^'^^ 

and sin a = Then, up to c^^ we have 

-1 dU -1 GM 1 x-x* 



yrr^Sa; ^TT^ R^ + ^(^+-)^ 1 - 6gg ^{x -x*Y + {y- rf 
-1 GM 



■ ^2 _|_ G{M-rym)R 



cosa, 



I.e. 



-1 GM 
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and analogously 



GM 



G(M-bm)R 
7^ 



az = Wx = Wy = Q 
GM i - 



sin a, 



GM 



Wz 



-I u' 



Gm VxD — VyX 



G{M-5m)R 



COS a = 



2 ^2 ^ G(M-5m)fl 



In the last equality it is used that cos a « ^, sina « ^, ^ and Uy « ^'y-^- Wc may assume 
approximately that x* /y* = x/y, because the barycenter is very close to coordinate origin. Hence 
cos a = x/r and sin a = y/r. 

By replacing these values in (2.18a,b), and neglecting the terms of order c~* and smaller, we obtain 



fx 



GM x — X* 5 GM xvx + yvy 2vy Gm xvy — yvx 

+ o ^ I {7.2a) 



R2 I G{M-5m)R R ' 2 B? 



c2 B? 



GM 



y — y* 5 Vy GM xVx + yvy 2vx Gm xVy — yVx 



R-i I G{M-^,n)R Jl ' 2c2 i?2 



C2 i?2 



, (7.26) 



These equations are analogous to (2.2), and continuing to deal like for the perihelion shift in section 
2, we obtain 

d / 2d'f\ cPy <Px _ 
_ 5 GM XVy - yvx xVx +yvy Gm xvx + yVy xVy - yVx _ 



2 i?2 



+ 2. 



(^M + 2m) 



G XVx + yVy XVy - yVx 
E? (? r 



B? -2 
5 



(^M + 2m 



r 

Gr"^ 



c2i?2 



V dt)dt\r) 



= -{^M + 2m 



GM^ 



2 dip \ dp 

' c^{M + mf \^ HJH' 
where p = 1/r. The solution of this differential equation is given by 



dt 



M + 2m 



G M2 

c2 (M + to)2 



(7.3) 



and hence 



dt ) 



G M"^ 

1 + (5M + 4m)^-— -^p 

c2 (M + m)2 ^. 



Analogously to (2.5), now we have 



/dp\ 
\dipJ 



92 = z;2C-2 



. . M^ 

1 + (5M + 4m)^-— -^p 

^ ^c2(M + m)2'^. 



(7.4) 



Further we should find v as a function of -R or r. We can replace 
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into the terms of and ay. According to this approximation, the solution of (2.18d) is given by 

\ c^) V (M + m)(M-5m) '^V M + m ) ~ 

+ 



( 



m^N-i Const. 



Hence after some transformations we can express v as a function of p up to c , and obtain 



^2_2 GM^p G'-M'p-(UI'- - AdiM - ()m^) 



(M + m)2 c2(M + m)4 

Note that the constant of integration Const, is neglected in the last formula, because one can verify 
that it has no role in the final formula for the periastron shift, like for the perihelion shift in section 
2. Substituting this value of into (7.4) we obtain 

{^y + p' = A + Bp + Cp\ 

where 

G^M^ G2m4(4M2 - 4mM - 6m^) _ 6G^M^ 



C = (lOM + 8rn) 



Cl{M + m)^c'^ Cl{M + mfc^ Cl{M + mfc'^' 



Since C2 = pTra^-^l — = Vl — e^, where is the semi-major axis of the relative orbit 

and P is the orbit period, we obtain 

3P^(M-Fm)^ 
27r2c2 af{l - e2) ' 

and analogously to the calculations for the perihelion shift in section 2, for the periastron shift in this 
case we obtain 

3G2(M + m)2p2 

Using that af = P'^G{M + m)/{4Tv^), we obtain finally 

67r(47r2)V3G2/3 

= p2/3,2(i _ ,2) + ™)'^'- (7-6) 

This value (7.6) of the periastron shift is the same which predicts the GR. 

Note that if we did not take into account that the pulsar with mass m has an influence on the 
gravitation from the companion with mass M, i.e. if we omit the expression 1 + in (7.1), then 
instead of (7.6) the periastron shift would be jj^^ At^, where Ay> is given by (7.6). It means that the 

companion star participates with j^^-pp^Aip in the total periastron shift given by (7.6). Symmetrically, 
the pulsar with mass m, participates with Aip in the total periastron shift. 
Now let us assume that 

-1 GM X , -1 GM V 

- — , and 



y'l _ J?2 I G(aM+0m)R r ' ^ ^/^ _ ^ j;2 , G(«M+/3m)ii ^ " 

Then the previous calculation leads that the periastron shift instead of (7.6) is equal to 

6.(4.2)1/3^2/3 (7 - a)M2 + {8- a- P)mM + (1 - /3)m2 
p2/3,2(i_,2) (M + m)/ X 6(M + m)2 " 
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Since the angle Aip should be symmetric with respect to M and m, it must be 7 — a = 1 — /3, i.e. 
a — (3 = But this is satisfied if and only if the "correcting multiplier" is 1 + ^r^. Further, in that 

case 

and it depends on the total mass M + m, but not on M and m separately. Finally, a must be 1 (and 
hence /? = — 5) in order to be 

6.(4.^)V3gV3 

^ p2/3c2(l _ e2-) 

for m = 0. 

Note that the formula (7.3) does not depend on a and /3. Thus, (7.3) for the pulsar's orbit can be 
important test for ILLTG, like the tests 1 and 2 in section 5. For example, from (7.3) it just follows 
that 

/2d^\ :(r^^) =1 {5M + Am)GM\ 

V dt J periastron V dt J apastron ttpC^ {M + m)^ {1 — c'^) ' 

where Gp is the semi-major axis of the pulsar's orbit. Using the recent data about the Hulse- Taylor 
pulsar PSR B1913+16 [12], wc get 



V dtJ periastron ' V dtJ 



= 1-0,47 X 10" 

apastron 



Probably using the fact that A(p should be a symmetric function with respect to M and m, or 
some similar conditions, it can be used to obtain the gravitational potential in the form 

GM C^MiM - 5m) _ G^MjaM'^ + (3mM + 7771^) 

and hence to obtain the post-Newtonian approximations of 2,3,... order. For post-Newtonian approx- 
imations of order 2,5; 3,5;... should be considered the multiple moments. 



8 ILLTG and the Hubble constant 

The representations in sections 1-7 would not be complete if we skip some cosmological consequences 
considered in [10,11]. In this section the basic arguments of [10] from viewpoint of ILLTG will be 

considered. 

In [10] is given a new view of the Hubble law, showing that the Hubble red shift and the anomalous 
acceleration of Pioneer spacecraft (see Ref. [1,2,3]) are consequences of the change of the potential in 

the universe. It is detected red shift such that = I'o^l — ^j, and replacing v = RH, indeed it is 

observed that 

u = uo{l-^). (8.1) 

On the other hand, if the potential changes linearly with the time, after time t = R/c we have 
additional potential ^ , and hence we have a frequency shift such that 

1 RdV\ 



u = 1/0(1 



(8.2) 



We should remark here that the sign of the potential V is such that near to the spherical gravitational 
body the potential {V = ^^) is larger than the potential far from gravitation (V = 0). 
Comparing the formulas (8.1) and (8.2) we obtain 
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The change of the potential energy in the universe is probably related with the dark energy, which is 
about 67% of the whole energy in the universe [6] . Some recent measurements show that H is about 
72 km/s/Mpc, but we shall accept value of about 82 km/s/Mpc. 

Since ^ < 0, in the past the time in the universe was slower than now, and in future it will be 
faster than now. Indeed, each second the time is faster approximately A = (1 + il x Is) = 1 + 3 76x10^^ 
times. 

On the other hand, in the last decades are considered the motions of the spacecraft Pioneer 10 and 

11, by comparing the initial frequency of a signal which is sent from the Earth to the spacecraft and 
the frequency of the signal which comes back. The frequency of the received signal on the Earth does 
not fit with the predicted frequency modeled by using the Doppler effect, the position and motion 
of the spacecraft using the GR and a lot of perturbations [3]. The observations show that appears 
almost a constant apparent unmodeled acceleration ap ~ (8, 74 ± 1,33) x 10~*cm/s^ which seems to 
act to the spacecraft toward the Sun. 

The acceleration ap is introduced by formula (15) in [3], i.e. 

[Vobs{t) - Vmodel{t)]DSN = -fQ— ^, whcrC Vmodel = Vq I - 2 ^"'°'^'''^^'^ , (8.4) 

c L c J 

where is the reference frequency, the factor 2 is because of two way data, vmodel is the modeled 
velocity of the spacecraft due to the gravitational and other forces, and i>obs is frequency of the re- 
transmitted signal observed by DSN antennae. After time 2t has been detected a small blue shift on 
the top of the red shift caused by the motion of the spacecraft outwards the Sun. Form (8.4), an 
unexplained blue shift 

. = .0(1 + ^), (8.5) 

is detected, where 2t is the time of the light signal in two directions. Without using the acceleration 
ap, the anomalous frequency shift can be interpreted by "clock acceleration" —at = —2, 8 x 10~^*s/s^ 
[3]. This causes just the blue shift given by (8.5). This model assumes that the time is uniform, but 
there is only technical problem with the clocks. Although this model gives good explanations for both 
frequency shift and the trajectories, it is rejected [3]. 

Many people have noticed that the acceleration ap and the Hubble constant H are related by [3] 

ap = cH. (8.6) 

Indeed, if we assume that iI=82km/s/Mpc, then from (8.6) it follows op = 8 x 10~^cm/s^. Thus we 
assumed that iJ=82km/s/Mpc. 

Now we shall explain the apparent acceleration ap connected with the frequency shift, according 
to ILLTG. We denote by x,y,z,t the orthonormal coordinates and denote by X, Y, Z, T ordinary 
coordinates where there is time dependent gravitational potential V = —(?Ht. According to (1.1), 
where V = —Hc^t, we have the following equalities 



dx 



= (1 + ^) ^dX = (1 + tH^X, dy=(l + ^) ^dY = (1 + tH^Y, 



dz=(l + ^) ^dZ = (1 + tH^Z, dt=(l- ^) ^dT = (l - tH^dT. (8.7) 

Rem,ark. Note that it is sufficient in this section to assume that (8.7) are satisfied. Then it is 
not necessary to speak about the time dependent gravitational potential. The coefficients 1 + Ht and 
1 — Htin (8.7) should be exponential functions, but neglecting H'^ and smaller quantities we accept 
these linear fimctions. 

If we use orthonormal coordinates x, y, z, t, then the motion of the spacecraft does not depend on 
the time dependent potential —Hc^t. Indeed, it is easy to see that the tensor Fij does not depend 
on this potential, and hence the motion of the particles (2.2) does not depend on the time dependent 
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potential, which means that ^) = (0, 0, 0), if we consider only the time dependent potential. 

Using the equalities (8.7), (2.2') and (8.6) we have an additional 3-acceleration 

'(PX (fY £Z\ ^/„dX „dY „dZ\ 



dTJ 



In the explanation of the shift frequency of Pioneers, the first component —3Ht(^^^, is 
smaller, compared with the second and main component 

/ „dY rrdZ\ /apdX ap dY apdZ\ 

\ dr'^dT'^dfJ^'^Tdr'Tdr'TdT)- 

We shall explain why it is measured the blue shift given by (8.5), instead of the red shift like in 
(8.1). It will be explained without using the acceleration (8.8), but only the definition (8.4) where ap 
is introduced, and our model of time dependent gravitational potential. 

Let R = R{T) he the distance from the spacecraft to the DSN antenna. Let us calculate ^ = 
(t'obs — ^o)/'^o- According to our model, neglecting the terms containing H^, this expression is equal 
to 

= -2H 2-^, 

va c c 

where the expression —2H^ corresponds to the Hubble red shift from (8.2). On the other side, 
according to (8.4), where Vmodei it determined without using the Hubble constant H, we have 

^ = 2«p4-2itl£^. 
vo c 

Prom the last two equalities we are able to find the expression for ap 

/IdR 1 /dR\ \ 

''^ = -iRdT-RydT)H=o'-'')- 

Using that ^ = ^(1 - 2HT) = (^) (1 - 2HT) according to (8.7), we obtain finally 

The equality (8.9) gives the required expression of the acceleration ap. Now we are not able to 
determine the initial value of T, i.e. when we should start to measure the time T in (8.9). It depends 
on the estimation of the initial values for the motion of the spacecraft, i.e. it depends on the departure 
of the real initial values. Thus, for different spacecraft are obtained close but different values of ap. 
Note also that for each initial value of T, when T tends to infinity (or R tends to infinity), then ap 
tends to cH. In an ideal case, if the experiment is done in an inertial system, then dR/dT = R/T is 
a constant velocity and then ap = cH. 

The change of the gravitational potential causes slight changes in the results of the known ex- 
periments about GR, but they arc negligible. The change of the results in the test of Shapiro time 
delay is negligible if we consider radio signals on short distances like in the solar system. Indeed, the 
change of the potential in the universe is about 2,4 x lO^cm^/s^ each second, while the gravitational 
potential on lAU from the Sun is equal to 30^km^/s^ = 9 x lO^^cm^/s^. This potential difference 
will be achieved from the universe after 118 years, which is too long compared with the time needed 
for an experiment about Shapiro time delay. The change of the deflection of the light rays near the 
Sun is completely negligible up to c^^. The angle between two pcrihclions is changed additionally of 
order ^ ^ , where Q is the orbit period of the planet and e is the eccentricity of the elliptic orbit. 
This change of the angle is negligible with respect to the measured angle ^^I^iL^i^ ■ 
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Note that if we replace V hy V ~ c^Ht {tH « 1) in the metric of GR 



= diag(^-l + _ - 1 + _, 1 + _, 1 + _ j 

the additional acceleration is not convenient. The situation becomes much worse if t is very large. 
While the Newton's theory is invariant under the gauge transformation V V + C, and ILLTG 
is invariant under the gauge transformation fj, i-^ C^, GR does not posses similar property. The 
equations of motion are unchanged if we replace the metric gij by Cgij, but it can not be correlated 
with the change of the gravitational potential. 
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